Coupling of spin and orbital excitations in the iron-based superconductor FeSe sTeo.s 
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We present a combined analysis of neutron scattering and photoemission measurements on su- 
perconducting FeSeo.5Teo.5- The low-energy magnetic excitations disperse only in the direction 
transverse to the characteristic wave vector (|,0, 0), whereas the electronic Fermi surface near 
(|,0,0) appears to consist of four incommensurate pockets. While the spin resonance occurs at an 
incommensurate wave vector compatible with nesting, neither spin-wave nor Fermi-surface-nesting 
models can describe the magnetic dispersion. We propose that a coupling of spin and orbital corre- 
lations is key to explaining this behavior. If correct, it follows that these nematic fluctuations are 
involved in the resonance and could be relevant to the pairing mechanism. 



The quest to understand the mechanism of high tem- 
perature superconductivity gained new momentum with 
the recent discovery of Fe-based superconductors Ex- 
perimental studies of these materials have shown that by 
chemically tuning the carrier density one can obtain both 
magnetically ordered and superconducting phases P The 
superconducting transition temperature is maximized for 
conditions close to where magnetic order is suppressed, 
leading to predictions that magnetic fluctuations are 
important for electron pairing and motivating compar- 
isons with copper-oxide superconductors. Initial theories 
have treated the magnetic excitations as independent of 
the atomic orbital character of the conduction electron 
statesj^an assumption that seems to work well in other 
metallic magnetic systems; however, there have been re- 
cent proposals that coupled spin and orbital order occur 
in the antiferromagnetic stateP^SD 

In this paper, we present a study of low-energy 
magnetic and electronic excitations in superconducting 
FeSeo.5Teo.5- Particularly striking is the anomalous 
anisotropic dispersion of magnetic excitations, which is 
quite distinct from the spin- wave- like excitations typi- 
cally seen in magnetic metals. Angle-resolved photoemis- 
sion spectroscopy (ARPES) reveals an unexpected split- 
ting of larger Fermi-surface pockets into smaller pockets 
with distinct orbital characters. A simple nesting picture 
does not appear to be compatible with the magnetic dis- 
persion. We propose, instead, that a strong coupling of 
spin and orbital correlations can lead to the observed 
anisotropy. Such a coupling would suggest that elec- 
tronic nematic fluctuations may be present in the nor- 
mal state and may participate in the superconductivity. 



This would be a rather different way to connect with the 
cuprates!^ 

Before continuing, let us first resolve our choice of co- 
ordinates. The structure of FeSeo.sTeo.s contains layers 
of Fe atoms forming a square lattice, with Se/Te atoms 
centered above or below these squares, alternating in 
a checkerboard fashion [Fig[TJa)]. The c-axis displace- 
ments of the Se/Te atoms break the translational sym- 
metry, doubling the unit cell size; however, as we will 
show, this effect is rather weak for the low-energy elec- 
tronic structure involving Fe-derived orbit als. Hence, we 
choose to work with a unit cell containing one Fe atom 
(a = b = 2.69 A and c = 6.27 A). Reciprocal lattice 
vectors will be specified in units of (27r/a, 27r/a, 2tt/c). 

It is already experimentally established that the 
low energy magnetic excitations in superconducting 
FeSe^Tei-^ occur near the in-plane wave vector Q = 
(|,0)P21 which is the same as the magnetic ordering wave 
vector found in .RFeAsO (R = La, Ce, Pr, Nd) and 
AFe2As2 (A = Ca, Sr, Ba)J^ where it corresponds to 
columnar antiferromagnetic order [Fig.[TJb)]. (Note that 
the magnetic ordering in Fei+^Te is different! 14 * 15 *) Sev- 
eral recent theoretical analyses have proposed that or- 
bital ordering, involving the Fe 3d xz and 3d yz states, is 
an essential feature of the magnetically-ordered statePtSU 
This is consistent with the fact that a symmetry-lowering 
structural transition, which breaks the degeneracy of 
the d xz and d yz states, always occurs before magnetic 
ordering.^ Very recent experiments emphasize the elec- 
tronic anisotro py of t his st at ep3IIHl demonstrating its ne- 
matic characte j 19 * 20 *. 

For the experimental studies, single crystals of 
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FIG. 1: (color online), (a) Atomic positions within an Fe 
plane of FeSeo.5Teo.5- Filled (open) diamonds indicate Se/Te 
atoms that sit above (below) the plane. Solid line: unit cell 
used here; dashed line: crystallographic unit cell, (b) Spin 
pattern corresponding to Q = (§,0). (c) Reciprocal space 
corresponding to (a), with symmetry point labeled. Solid line: 
Brillouin zone for one Fe per unit cell; dashed line: Brillouin 
zone for two Fe per unit cell, (d) Reciprocal space labeled in 
reciprocal lattice units, indicating directions of scans for the 
neutron scattering measurements. 



FeSeo.sTeo.s were grown at Brookhaven by a unidirec- 
tional solidification method. The nominal composition 
had no excess Fe, and T c , measured by magnetic sus- 
ceptibility, is 14 K. 21 The ARPES spectra were recorded 
on beamline U13UB at the National Synchrotron Light 
Sounce using a Scienta SES2002 electron spectrometer. 
The incident photon energy was 17 eV, while the photo- 
electron energy and angular resolution was 15 meV and 
0.1°, respectively. Samples (small pieces from the crystal 
studied with neutrons) were cleaved in situ (base pressure 
of 5 x KT 11 Torr) and held at T = 15 K for the mea- 
surements. Complementary ARPES data were collected 
on beamline 12.0.1 at the Advanced Light Source (ALS), 
using 50-eV photons (energy resolution of 25 meV) and 
T = 15K. 

The neutron scattering experiments were performed at 
the cold-neutron triple-axis spectrometer SPINS located 
at the National Institute of Standards and Technology 
Center for Neutron Research (NCNR). The 9-g crystal 
was aligned in the (hkO) scattering plane; a mosaic width 
of < 0.6° (limited by resolution) was measured. The 
incident neutron energy was selected with a vertically- 
focusing pyrolytic graphite (PG) (002) monochromator. 
Scattered neutrons were analyzed with five 2.1 cm x 15 
cm PG(002) flat analyzer blades that reflected neutrons 
with Ef = 5.0 meV onto a 3 He proportional counter. 
A cooled Be filter was placed between the sample and 
analyzer to suppress higher order neutron contamination. 

Let us consider the ARPES results first. The electronic 
spectral weight within 0-10 meV of the Fermi energy, Ef, 
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FIG. 2: (a) ARPES intensity map at Ef, integrated over 
0-10 meV binding energy; upper half measured at the ALS 
with 50-eV photons, lower half from the NSLS with 17-eV 
photons. Thick red lines indicate estimated Fermi-surface po- 
sitions; thin red lines denote pockets inferred by symmetry. 
Red arrows denote polarization direction; at (0, 0.5) polariza- 
tion has a k z component that varies with k x . (b) Intensity 
map (white — minimum, blue — maximum) along line de- 
noted as Cut A in (a). Red dots indicate positions of bands 
as they cross Ef- Asymmetry about k y = is due to imper- 
fect sample orientation, (c) Calculated spectral weight at Ef, 
averaged over k z . Orbital coding: red-c^; b\ue-d xz ; green- 
all others, (d) Calculated band dispersions along symmetry 
directions, with intensity scaled according to spectral weight, 
obtained by unfolding bands calculated in a reduced zone^ 
color code same as (c). 



is plotted in Fig. [5|a); the positions of identified band 
crossings, determined from energy- momentum cuts as in 
Fig. [2jb), are indicated by thick red lines. For compar- 
ison, our local density calculation of the Fermi surface 
is shown in Fig. 2(c) and the back- folded band struc- 
ture in Fig. 2(d). Several features are obvious. First of 
all, we can see that the effect of the unit cell doubling 
due to the Se/Te positions is small. To appreciate this, 
note from Fig. [TJc) that for the doubled unit cell the T 
and M points correspond to the same crystal momen- 
tum. While the band dispersions must be the same at 
both points, the spectral weights of the individual bands 
need not bep^and that is certainly the case here. Indeed, 
in Fig. |2|d) one can barely recognize the relationship be- 
tween M and T points, reflecting weak Umklapp effects 
on the Fe-orbitals. 

Secondly, experiment indicates the presence of small, 
incommensurate pockets about the X and Y points, in 
contrast to the larger commensurate pockets predicted by 
our calculation. The spectral intensities of these pockets 
are quite sensitive to the orientation of the photon po- 
larization. We can understand this dependence as fol- 
lows. The theoretical pockets about X (and Y) have 
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distinct orbital character, as indicated in Fig. [2|c) and 
(d). Fe d yz character (red) dominates near (|,±fco,0), 
d xz (blue) near (=b/io, ^ , 0) , and d xy around (| — /io>0, 0) 
and (0, - — fc,0). Due to the broken cubic symmetry at 
the Fe sites and the hybridization with the neighboring 
Te/Se p orbitals, d yz /d xz /d xy Wannier orbitals actually 
contain p x /Py/Pz character (see Ref. [6]for a related anal- 
ysis of Wannier functions). Assuming that the ARPES 
matrix elements are primarily p — >• s (which is reasonable 
for low photon energies), the association between orbital 
character and polarization direction is: d yz : [100], d xz : 
[010], d xy : [001]. From the observed intensities, we con- 
clude that the pockets at (|, ±fco, 0) are d^-like, those at 
(=b/i , § 5 0) are d^-like, and those at (0, \ ± fco>0) [and, 
by symmetry, at (| ± /io>0,0)] are d^-like. The orbital 
character of the spectral weight at Ep is qualitatively 
compatible with calculations P 2 ^ At the same time, the 
observed small hole pockets are different from the larger 
Fermi surfaces of the local density calculation, probably 
due to the inability of this mean-field technique to prop- 
erly account for dynamical magnetic correlations. 

Now let us turn to the neutron scattering results. As 
already mentioned, a spin resonance has been observed 
to develop below T c in FeSe^Tei-^ with x = 0.4 at 
the in-plane wave vector Q = (|,0) and an energy of 
6.5 meVP^ We have confirmed the presence of a reso- 
nance in our x = 0.5 sample. What we wish to focus on 
here is the Q dependence of the magnetic response in the 
vicinity of the resonance. In the superconducting phase 
(T = 1.5 K << T c ), transverse scans along (0.5, k) [scan 
A m Fig. [ljd)] exhibit pairs of peaks at finite \k\ that 
disperse with energy above 5 meV, as shown in Fig.|3ja). 
The color-coded map of intensity vs. Q at the resonance 
energy of 6.5 meV, Fig. |3jb), demonstrates an intrigu- 
ing anisotropy: the transverse peaks are not reproduced 
along the longitudinal (0.5 + /i, 0) direction as one would 
expect from Fermi surface nesting between a pocket at V 
and the four pockets about X. 

In the normal state (T = 20 K), the spectral weight 
of the resonance is moved to lower energies, and the 
peaks appear to remain split down to 4 meV [blue sym- 
bols in Fig. |3ja)]; the asymmetry between transverse 
and longitudinal directions remains the same as below 
T c . The dispersion of the peaks along the transverse di- 
rection is plotted in Fig. |3jc), with points at 10 meV 
and above from Lumsden et al. 2A (For similar recent re- 
sults, see Refs. 25 26.) The line through the points is 
a fit to hjj = Eos'm(7r(k — ko)) with Eo = 121 meV for 
k < ko = —0.09 . If th ese excitations were like spin waves, 
as in CaFe2As2 p we would expect to see cone-shaped 
dispersions coming out of (0.5, ±&o), which is clearly not 
the case. The dispersion of isolated intensity peaks along 
a single direction is quite unusual and requires considera- 
tion of factors beyond the degree of electronic correlation. 

We propose that the solution to this puzzle involves 
the coupling of spin and orbital correlations. As previ- 
ously mentioned, such a coupling has already been pro- 




FIG. 3: (color online). Momentum dependence of the spin 
resonance and low energy spin fluctuations, (a) Inelastic neu- 
tron scattering intensity as a function of Q, obtained with 
constant energies, 8 meV, 6.5 meV, 5 meV, and 4 meV, at 1.5 
K (< T c ) (open red circles) and 20 K (> T c ) (filled blue cir- 
cles), taken along the transverse (0.5, k) direction. The lines 
are fits to Gaussians; diamonds (open, 1.5 K; filled, 20 K) indi- 
cate fitted peak positions. Horizontal bars represent the full- 
width-of-the-half-maximum of the instrumental Q-resolution. 
(b) Color contour map summarizing 1.5-K data obtained from 
a series of (0.5, k) scans centered at positions marked by red 
squares, (c) Dispersion of the spin fluctuations. Open green 
squares extending to high energies are taken from Ref. 
while open red diamonds (1.5 K) and filled blue diamonds 
(20 K0 at low energies are our data taken from (a); hori- 
zontal bars are fitted peak widths at the lowest energies and 
T = 20 K. The line is explained in the text. 



posed for the antiferromagnetic phaseP^In fact, it has 
been showrP that the reduction of lattice symmetry can 
be associated with long-range orbital ordering, and the 
resulting anisotropy of the superexchange paramet ers is 
consistent with the observed spin-wave dispersions! 27 * 28 * 
In the present case, the orbital degeneracy of d xz and 
d yz orbitals re-activates the orbital freedom in the ab- 
sence of long-range orbital order. This changes entirely 
the nature of the low-energy spin excitations. 

Spin correlations characterized by Q ~ (^0), as 
observed here, correspond to the columnar pattern of 
Fig. |TJb) . As shown in Ref. |6j such correlations are di- 
rectly tied to the spin polarization of the d yz orbital and 
its ferro-orbital correlation. [Similarly, spin correlations 
at Q « (0, \) are tied to the d xz orbital.] We propose 
that the spin-flip excitations in the present case are con- 
strained by a strong coupling to orbital correlations, with 
two leading contributions illustrated in Fig. [4j First, if 
we consider flipping a single spin, we find that the or- 
bital also needs to flip from d yz to d xz [c.f. Fig. j4^a)], in 
order to better utilize the superexchange energy. (If the 
orbital remains unchanged, a local spin flip would have 
closed almost all the superexchange paths within the d yz 
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FIG. 4: (color online), (a) Cartoon of spin+orbital flip excita- 
tion discussed in text, assuming finite range spin and orbital 
correlations. Vertical (horizontal) line denotes d yz (d xz ) or- 
bital, (b) Form factor for the excitation in (a). Color map: 
low-medium-high = blue-green-red (dark-light-gray), (c) 
Cartoon of inter-site spin-flip excitation, (d) Form factor map 
for the excitation in (c). 



and d xz subspace.) This "inter-orbital spin/orbital flip" 
process introduces a strongly anisotropic form factor that 
goes to zero along the longitudinal direction, as shown in 
Fig. [4^b). This is consistent with the strong anisotropy 
of our neutron data and the previous observation at high 
energies p2 Another relevant process, one associated with 
itinerancy, is the "inter-site spin flip" process shown in 
Fig. N|c); the latter has a large form factor around (|,0) 
[c.f. Fig.Qd)]. Interestingly, due to the columnar AF ar- 
rangement of the short-range spin correlations together 
with the Pauli exclusion principle, our proposed inter- 
site spin-flip process can only propagate along the y- 



direction at low-energy, in good agreement with the ob- 
served anomalous dispersion. A quantitative analysis of 
the dispersion (beyond the scope of this paper), will re- 
quire evaluating the excitation energy as a function of 
wave vector, taking into account both processes. In any 
case, the essence of our proposal is that the spin excita- 
tions are hybrids of the magnons and orbitons that have 
been postulated previously!^ 

In conclusion, we have presented experimental evi- 
dence that the electronic and magnetic excitations of 
superconducting FeSeo.sTeo.s are different from those of 
common band- structure and spin- wave models. We have 
proposed that the intensity anisotropy of the magnetic 
excitations may be explained by considering coupled spin 
and orbital excitations. If correct, this would imply that 
nematic excitations are a key feature of the normal state 
from which the superconductivity develops. While the 
details of the magnetism are quite different, the nematic 
correlations would provide an intriguing connection to 
the physics of cuprate superconductors.^ 
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